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m ■ Abstract 
> ■ 

m 

■ We study constraints on the CP violating phase 7 in the Kobayashi-Maskawa 



model using available experimental data. We first follow the conventional 
method to up date the constraint on 7 by performing a analysis using data 
from {exl, ^"n^Bas \^ub/ycb\- We also include the recent information on 



Oh: 

^ • sin 2/? in the analysis. We obtain the best fit for 7 to be 66° and the 95% 



C.L. allowed range to be 42° ~ 87°. We then develop a method to carry 
out a analysis based on SU(3) symmetry using data from B vrvr and 
B — > Kir. We also discuss SU(3) breaking effects from model estimate. We 
find that present data on i? — > vrvr, Kir can also give some constraint on 7 
although weaker than the earlier method limited by the present experimental 
errors. Future improved data will provide more stringent constraint. Finally 
we perform a combined fit using data from {exl, ^n^Bds^ l^ub/Vcbl, sin2/3 
and rare charmless hadronic B decays. The combined analysis gives 7 = 67° 
for the best fit value and 43° ~ 87° as the 95% C.L. allowed range. Several 
comments on other methods to determine 7 based on SU(3) symmetry are 
also provided. 
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I. INTRODUCTION 



The origin of CP violation is still a mystery although it has been observed in neutral 
kaon mixing more than 35 years. One of the most promising model for CP violation is the 
model proposed by Kobayashi and Maskawa in 1973 ||l|]. This is now referred as the Standard 
Model (SM) for CP violation. In this model CP violation results from a non-removable phase 
7 in the charged current mixing matrix, the Cabbibo-Kobayashi-Maskawa (CKM) matrix 
|Jl|J^, VcKM- There are also other mechanisms for CP violation. To understand the origin 
of CP violation, it is important to study in every detail of a particular mechanism against 
experimental data. In this paper we carry out a study to constrain the CP violating phase 
in the SM using available experimental data. 

The CKM matrix Vckm is a 3 x 3 unitary matrix and is usually written as 



^ Vud Vus Vub ^ 



CKM 



V 



(1) 



Vcd Vcs Vcb 

Vtd Vts Vtk 

In the literature there are several ways to parameterize the CKM matrix. The standard 
particle data group parameterization is given by 0] 



CKM 



C12C13 



S12C13 



-S12C23 - Ci2S23Sl3e*'^" C12C23 - Sl2S23Sl3e*''^=* 



V 



■512-523 ~ Ci2C23Si3e 



"C12S23 — Sl2C23'Sl3e 



S23C13 
C23C13 



(2) 



where Sij = sin 6ij and Cij = cos dij are the rotation angles. A non-zero value for sin5i3 

violates CP. Another commonly used parameterization is the Wolfenstein parameterization 

^ which expands the CKM matrix in terms of A = \Vus\ and is given by 

\ 



/ 



CKM 



1 \2 



C(A^). 



(3) 



A A\^{p-ivi) 
-X 1 - iA^ AX' 

^AX^il- p-i7]) -AX' 1 
The parameters A, p, rj are of order unity. When discussing CP violation in kaon system, it is 
necessary to keep higher order terms in A, namely, adding —A'X^{p + iri) and —AX'^lp + irj) 



J 



to Vcd and Vts, respectively. CP violation in this parameterization is characterized by a 
non-zero value for rj. 

Due to the unitarity condition, one has 

v:,v^d + KtKd + V*Vu = 0. (4) 

In the complex plane the above equation defines a triangle with angles a = 
-Arg{VtdV;jVudV:,), (3 = -Arg{V,aV:JVuVg) and 7 = -Arg{VudV:JV,dV:,) as shown 
in Figure 1. 




V V / \ V V 

* ud * ub/ \* td * tb 
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FIG. 1. The KM unitarity triangle. 

To a very good approximation the phase (5i3 is equal to 7. In terms of p and 17, the angles 

a, (3 and 7 are given by: 

. 27]{p^-p + if) . 2?7(1 -p) _^r} 
sm2a = -- — — — — — -, sm2/3 = - — — 7 = tan -. 5 

[[I - PY + ri^)[p^ + rj^) [l-pY + rj^ p 

In this paper we will concentrate on obtaining constraint on the phase 7. Great efforts 
have been made to constrain or to determine the CP violating phase 7. Previous studies 
mainly used experimental data on: i) The CP violating parameter in the mixing of 
neutral kaons; ii) The mixing parameters Am^^ and Am^^ in Bd^s — B^^s systems; And 
iii) |\4fe/V^b| which characterizes the strength of the charmless fiavor changing and charmed 
fiavor changing semi-leptonic B decays. The best fit value for 7 from these considerations 
is around 65° 



During the last few years, several rare charmless hadronic B decays have been measured 
0. Some of these decays are sensitive to 7 and therefore can be used to constrain it P,^. 
Analysis based on naive factorization approximation suggests that 7 tends to be larger than 
90° in conflict with the analysis mentioned earlier ^j^. If confirmed, it is an indication of 
new physics beyond the SM. Of course due to uncertainties in the experimental data and 
theoretical calculations, it is not possible to draw a firm conclusion whether this conflict 
is real at present. To improve the situation, in this paper we will carry out an analysis 
replacing the naive factorization assumption by more general SU(3) flavor symmetry for the 
rare charmless hadronic decays of B to two SU(3) octet pseudoscalars Pi and P2, that is, 
B — > PP decays. 

SU(3) analyses for B decays have been studied by many groups and several interesting 
results, such as relations between different decay branching ratios, and ways to constrain 
and/or to determine the phase 7, have been obtained [p!0|-p^. SU(3) symmetry is expected 
to be a good approximation for B decays. At present experimental data from B —>■ DK{tt) 
support such an expectation |T^. However more tests are needed, especially in rare charmless 
hadronic B decays. Recently it has been shown that such tests can indeed be carried out for 
rare charmless hadronic B decays in an electroweak model independent way in the future 
|13| . Before this can be done, however, SU(3) symmetry can only be taken as a working 



hypothesis. In the rest of the paper we will study constraints which can be obtained from 
rare charmless hadronic B PP decays based on SU(3) symmetry. We will also study 
SU(3) breaking effects using model calculations. 

The paper is arranged as follow. In section II, we will review and up date the constraint 
on 7 using information from ex, ^f^B^s ^"^^ iKib/Kbl, and also information from sin 2/5 
measurement. In section III we will carry out a analysis of 7 using rare charmless hadronic 
B PP decay data based on SU(3) symmetry. We will also discuss SU(3) breaking effects. 
In section IV, we will make a combined study using results from sections II and III. And 
in section V, we will discuss some of the implications of the results obtained and draw our 
conclusions. 



II. CONSTRAINT ON 7 FROM |eA-|, AMb^^s, \Vub/Vcb\ AND sin 2/? 

In this section we first review and up date the constraint on 7 using experimental and 
theoretical information on e^:, Am^^^ and |V^;,/V"cb|. Such an analysis has been carried 
out before. The analysis in this section is an up date of the previous analyses which also 
serves to set up our notations for later use. We then include experimental data from sin 2(3 
measurement into the analysis to obtain the best fit value and allowed range for 7. 

There exist quite a lot of information about the CKM matrix |^. The value of Vus is 
known from Ki-^ decay and hyperon decays with good precision: 

A = 0.2196 ±0.0023. 

The parameter A depends on A and on the CKM matrix element \Vch\- Using experimen- 
tal data from B D*l^v and B —>■ Dl^u and inclusive b cW, analysis from LEP data 
obtains Kfe = 0.0402 ± 0.0019, and data from CLEO obtains = 0.0404 ± 0.0034. The 
central values of these two measurements are close to each other. In our analysis we will use 
the averaged value which leads to A = 0.835 ± 0.034. 

The value for \Vub\ has also been studied using data from B vr/z/^, B — > pWi and 
inclusive b —>■ uWi with 

\Vub/V,b\ = X^Jp^ + = 0.090 ± 0.025 (6) 

To separately determine p and 77 (or 7), one has to use information from other data. In the 
rest of this section we will carry out a analysis using constraints from the measurements of 
{exl, ^f^Bi , and iKife/Kbl along with other known experimental and theoretical information. 

The parameter ex indicates CP violation in neutral kaon mixing. The short and long 
lived mass eigenstates Ks and Kl of the neutral kaons can be expressed as the linear 
combination of weak interaction eigenstates K*^ and as \Ks) = p\K^) + q\K^) and 
\Kl) = p\K^) — q\K^). p and q are related to the CP violating parameter ex in Kl decays 
by: 
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q I -ex 

The precise measurements of the Ks — > tttt and Kl — > tttt decay rates imply ^ : 

leA'l = (2.271 ±0.017) x 10^1 
Evaluating the so called "Box" diagram, one obtains 

\eK\ = ^fj^^^f^Bj,iA'\^v) [y.ivctfsiyc, Vt) - Vcc) + Vuytf2iyt)A'\\l - P)] . (8) 



where T]tt = 0.574±0.004, rj^ = 0.47±0.04 and rj^c = 1.38±0.53 are the QCD correction 



factors, ArriK = rriK^ - uiks = (0-5300 ± 0.0012) x 10-2ps"\ and Bk = 0.94 ± 0.15 ^ is 



the bag factor. The functions /2 and /s of the variables yt = ml/m^ and yc = ml/m^r are 
given by p!7| : 

1 9 3 3x^\nx 



f2{x) 



4 4(1 -x) 2(1 -x)2 2(1 -x)= 



Neutral mesons and B^ show a behavior similar to neutral kaons. The heavy and 
light mass eigenstates, Bl and Bh, are different from B^ and B^ and are given by 

\Bl) = p\B',) + qm , 

\BH)=p\B',)~qm. (10) 

The mass difference Am = fTLB^ — friBj^ can be measured by means of the study of the 
oscillations of one CP eigenstate into the other. The world average value for Auib^ is |jl8 



Ams, = 0.487 ±0.014 ps-\ (11) 

The contribution to AfnB^ is from analogous "Box" diagrams as that for ex, but with 
the dominant contribution from the top quark in the loop. One obtains 



G 



2 



A^B, = ^ml,mBAfB,^BB,rVBytf2{yt)A'X'[{l - pf ± r,\ (12) 



where fs^^BB^ = 0.215 ± 0.040GeV r]B = 0.55 ± 0.01 []T5[ and the function /a is given 
by Eq. i). 

and mesons are beheved to undergo a mixing analogous to the B^ and B^. Their 
larger mass difference Aitlb^ is responsible for oscillations that are faster than the B^ and 
B^ oscillation, and have thus still eluded direct observation. A lower limit has been set by 



the LEP, SLD and CDF collaborations, as |18 



as: 



Atub, > 14.9ps~^(95% C.L.). (13) 
The expression for Auib^ in the SM is similar to that for ArriB^- Aitlb^ can be written 

AruB. = Am^.^^e^- (14) 

where all the theoretical uncertainties are included in a quantity ^, which is given by ||19|: 



e= f^'^^^ = 1.14 ±0.06. (15) 

fBa\ BBa 



The p and rj parameters can be determined from a fit to the experimental values of the 
observables described in the above. In the analysis we will adopt the strategies used in 
previous analysis in the literature fixing the known parameters, theoretical or experimental, 
to their central values if their errors were reasonably small reported in the left half of Table 
I. The quantities affected by large errors will be used as additional parameters of the fit, but 
including a constraint on their value as shown by the right half of Table I. All errors will 
be assumed to be Gaussian. This assumption may result in stringent constraints more than 
actually can be achieved because some of the errors may obey different distributions, for ex- 
ample those errors come from theoretical estimates may obey fiat distribution. Nevertheless, 
the results provide a good indication for the values of the parameters involved. 

To obtain the best fit values and certain confidence level allowed ranges for the relevant 
parameters, we perform a analysis using the above information. The procedure for 
analysis here is to minimize the following expression: 



2 _ (A - Af {rn, - {fut - rritf {Bk - BkY {fhc - Vcc? 



a, 



mt 



a 



Bk 



+ 



+ 



a. 



+ 



+ 



Vet 



+ 



Vcc 

■ \V^\ _ \Vub\\2 
• \V,b\ \Va\ > 



a 



\Vub\ 



+ 



cr 



+ x'(A(AmBj,aA(AmBj). 



AmB, 



(16) 



The symbols with a hat represent the reference values measured or calculated for given 
physical quantities, as listed in Table I, while the corresponding a are their errors. The 



parameters of the fit are p , t], A, rric, rrit, Bk, r]ct, t]cc, fsax Bb^ and ^. 



TABLE I. Input parameters for analysis using data from ex, ^^Bds ™^ iKtfe/K 



cb\ 



Fixed values 



Varied parameters 



A = 0.2196 ±0.0023 

Gp = (1.16639 ± 0.00001) x 10"^ GeV^^ 

fx = 0.1598 ibO.OOlSGey 

AruK = (0.5300 ± 0.0012) x 10"^ pg-i 

rriK = 0.497672 ± 0.000031 GeV 

mw = 80.419 ± 0.056 GeV 

mBa = 5.2794 ± 0.0005 GeV 

mB, = 5.3696 ± 0.0024 GeV 

ryj4 = 0.574 ± 0.004 

7/B = 0.55 ±0.01 
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A = 0.835 ± 0.034 

ryct = 0.47 ± 0.04 

rice = 1.38 ± 0.53 

mc(mc) = 1.25 ±0.10 GeV 

mt{mt) = 165.0 ± 5.0 GeV 



fs.VB^, = 0.215 ± 0.040 GeV 

Bk = 0.94 ±0.15 

^ = 1.14 ±0.06 

\€k\ = (2.271 ±0.017) X 10-3 

AruBd = 0.487 ±0.014 ps^^ 

\Vub/Vcb\ = 0.090 ±0.025 



i 

il 

i 

i 

m 
m 
m 
i 

El 



The inclusion of the ArriB^ data needs some explanation. The experimental data consists 
of measured values of A^AtjibJ and cr^(AmBj for various values of AniB^ plot in Figure 
2. To include this data in the fit, for each set of free parameters {A, p,ri,^) we calculate the 
value of ArriB, and find the corresponding experimental values of A and cr^ in Figure 2. A 
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nonzero value of Atub^ implies that there is — mixing, if observed one should have 
^ = 1 and otherwise ^ = 00. We follow Ref. to add to the total in Eq. (16) a 
Ax^ for the corresponding set of {A,p,ri,^), 

Ax' = x\AiAmsJ,aAAmnJ) = {^—^f. (17) 

Exp[—Ax''/2] is an indication of how likely a mixing with a given Aitib^ was measured by 
experiment. The sign of the deviation A — 1 should also be carefully treated. Naively the 
expression of Ax^ implies that a lower probability is attirubted to the Ams, values with 
^ > 1 with respect to Aitlb^ values having A = 1. To avoid this undesired behavior, we 
follow Ref. [^] to set A to unity for the range with A larger than one in Fig. 2. 



■a 
s 



< 



World average (prel.) 



data ± 1 a A 95% CL limit 14.9 ps 

1.645 a -I- sensitivity 17.9 ps"^ 

H data ± 1.645 <; 
I I data ± 1.645 a (stat only) 




2.5 



Am (ps ) 



FIG. 2. Experimental data on AmB^[17]. 

After p and 77 are determined it is easy to obtain the values of the angles in the unitarity 
triangle using relations in Eq. (|^) . The best fit values and the allowed regions in the p — r] 
plane are shown in Figure 3. The best fit values and their 68% CL. errors are 



p = 0.18^°:;^^, 77 = 0.34; 



-0.07 
-0.06 ; 



sin 2a 



-O.I9H2', sin2/3 = 0.70l°:J^, 7 



62+^2° 

"^-13° 



(18) 
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The 95% C.L. allowed regions for the above quantities are expressed as: 



0.03 < p< 0.38, 0.23 < < 0.50, 

-0.85 < sin 2a < 0.42 , 0.49 < sin 2/3 < 0.94 , 39° < 7 < 84°. 



(19) 



These results agree with previous analyses 
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FIG. 3. Constraints on p and r] using data from |e/<|, Am^^^ and iKib/Vctl, and sin2/3. In the 
figure on the left, only \eK\^ ^mB^s iKib/Vctl are used. The best fit value is indicated by the 
"+" symbol. The region in the dotted curve corresponds to ~ Xmin ~ f allowed region which is 
at the 39% C.L.. The 68% C.L. allowed region is within the dashed curve and the 95% C.L. allowed 
region is within the solid curve. The straight ray lines are the results for direct measurement of 
sin 2/3. The thick solid lines are for the central value of sin 2/3. There are two allowed regions. The 
region outside the two thin solid straight lines for each allowed region are excluded by the sin 2/3 
measurement at 95% C.L.. The 68% C.L. allowed regions are between the dashed lines. The figure 
on the right is a fit with sin 2/3 data also included in the ■ 

The solid line in Figure 4 is a plot of the minimal a function of 7 for the fit using 

1 6^:1, ArriB^^ and |V^tif,/Vcf,|. It is clear that changes with 7 dramatically. When going 
away from the minimal, raises rapidly indicating a good determination of 7. 

There are also direct measurements of sin 2/3 by several groups from the time dependent 
CP asymmetry in i? ^ J/ipKs- In the SM this asymmetry is given by 

V{B\t) ^ JlijKs) - nB\t) ^ J/,pKs) 



a{t) 



sin 2/3 sin( Am B^t) . 



(20) 
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FIG. 4. The solide line is the as a function of 7 using data from |e/^|, Am^^^ and 
The dashed Une is the as a function of 7 with sin 2/3 data included in the fit. 

The values measured by different groups are, 

0.34 ±0.20 ±0.05; 



sin 2(3 



n p:o+0.32+0.09. 
U-JO-0.34-0.10) 



BaBar [21] 
Belle P 



(21) 



CDF 22 



0.841?:^^ ±0.16 ALEPH 

The averaged value is sin 2/3 = 0.46 ± 0.16. 

For a given sin 2/5 there are, in general, four solutions for 7 with two of them having 
negative t] and another two having positive rj. To determine which one of them is the right 
solution, one has to use other information. Using the information from our previous fit, we 
can rule out some of the solutions. The allowed ranges for p and 77 from the averaged value 
for sin 2(3 is shown in the figure on the left in Fig. 3 by the straight ray lines. Since the fit 
from \eK\, ^f^B^s ^"^^ iKi^/Kbl determines > 0, only solutions with > are allowed. 
It is clear that one of the values for p and 7] determined from sin 2/3 measurement can be 
consistent with the fitting results in Eqs. ( pTSD and (p^). 

It is intersting to note that sin 2(3 data can eliminate a large allowed range in the p vs. r] 
plane at the 95% CL. level. One can also include the measured sin 2/5 into the analysis. 
The results are shown in the figure on the right in Fig.3. The as a function of 7 is shown 
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by the dashed hne in Fig. 4. The best fit values and their 68% CL. errors are given by: 

._nio+0.10 Q oQ+0.05 

sin2a = -0.19l|]i^, sin 2/? = 0.6ll[J:|J?, 7 = 661^^0. (22) 
And the 95% CL. allowed regions for the above quantities are: 
0.01 < p< 0.30, 0.21 <r]< 0.41, 

-0.88 < sin 2a < 0.45 , 0.40 < sin 2/3 < 0.80 , 42° < 7 < 87°. (23) 

III. DETERMINATION OF 7 FROM CHARMLESS HADRONIC B DECAYS 

In this section we study how the phase 7 can be constrained from experimental data on 
B PP decays, based on flavor SU(3) symmetry consideration. 

A. The Quark Level Effective Hamiltonian 

The quark level effective Hamiltonian up to one loop level in electroweak interaction for 
charmless hadronic B decays, including QCD corrections to the matrix elements, can be 
written as 

H!ff = %[K5k;(ciOi + C2O2) - E(K6Kvr + v,,v:^cT)o,]. (24) 

The coefficients ci,2 and c^'^ = — , with j indicates the internal quark, are the Wilson 
Coefficients (WC). These WC's have been evaluated by several groups ||2^, with |ci,2| >> 
|c^|. In the above the factor VcbV*^ has been eliminated using the unitarity property of the 
CKM matrix. The operators Oj are defined as ||24|| , 

Oi = {qiUj)v-A{ujbi)v-A , O2 = {qu)v^A{ub)v-A , 

03,5 = {qb)v-AEg'{q'q')vTA , 04,6 = {qibj)v-AEq'{qjq'i}vTA , 

(25) 

07,9 = |(<?&)v-AEq' eq>{q'q')v±A , 08,io = |(g»fei)y-A Eg' eq^q'^ql) 
Oii = ^ga^.G^-(l + 75)&. 
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)V±A 



where {qiq2)v-A = ^17^(1 ~ 75)'?2, G^'^ is the field strengths of the gluon, respectively. We 
have neglected the photonic dipole penguin term whose contribution to hadronic charmless 
B decays is small. The usual tree-level W-exchange contribution in the effective Hamiltonian 
corresponds to O2. Oi emerges due to the QCD corrections. The operators 03,4,5,6 are from 
the QCD-penguin diagrams. The operators O7,..., Oiq arise from the electroweak-penguin 
diagrams. On is the gluonic dipole penguin operator. 

The WC's at /i = 5 GeV with as{mz) = 0.118, in the regularization independent scheme 



in Ref . 



are 



where 



Cl = 


-0.313 , 


C2 = 


1.150 , 




0.017 , 


C4 — 


-0.037 , 




0.010 , 


r* - 


-0.046 , 




-O.OOlaem , 


r* - 






-1.321aem , 


r* - 


-- 0.267aem 




= -0.143 . 






are 


given in Ref. 
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(26) 



B. SU(3) Structure of the Effective Hamiltonian 

To obtain B decay amplitudes, one has to calculate hadronic matrix elements from 
quark operators. At present there is no reliable methods to calculate these matrix elements 
although simple factorization calculations provide some reasonable results for some decays. 



but not all of them [Q. It motivates us to carry out model independent analysis by studying 
properties of the effective Hamiltonian under SU(3) flavor symmetry and use them to obtain 
information about related decays. 

In general the decay amplitudes for B PP can be written as 



A{B ^ PP) =< P P\Htff\B >= ^[V^,v:^T{q) + V,,V*P{q)] 



(27) 



where T{q) contains contributions from the tree operators Oi,2 as well as penguin operators 
O3-11 due to charm and up quark loop corrections to the matrix elements, while P{q) 
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contains contributions purely from penguin due to top and charm quarks in loops. The 
amplitude T in Eqs. ( |27| ) is usually called the "tree" amplitude which will also be referred 
to later on in the paper. One should, however, keep in mind that it contains the usual tree 
current-current contributions proportional to Ci,2 and also the u and c penguin contributions 
proportional to c^'^ with 2 = 3 — 11. Also, in general, it contains long distance contributions 
corresponding to internal u and c generated intermediate hadron states. In our later analysis, 
we do not distinguish between the tree and the penguin contributions in the amplitude T. 

The relative strength of the amplitudes T and P is predominantly determined by their 
corresponding WC's in the effective Hamiltonian. For AS = charmless decays, the dom- 
inant contributions are due to the tree operators O12 and the penguin operators are sup- 
pressed by smaller WC's. Whereas for AS = —1 decays, because the penguin contributions 
are enhanced by a factor of VtbV^*/VubV*g ~ 50 |0| compared with the tree contributions, 
penguin effects dominate the decay amplitudes. In this case the electroweak penguins can 
also play a very important role [^. 

The operators Oi^2, 03_6,ii, and Oj^io transform under SU(3) symmetry as 3+3' + 6 + 15, 
3, and 3 + 3' + 6 + 15, respectively. We now give details for the decomposition under SU(3) 
for some operators. For AS = decays, O2 can be written, omitting the Lorentz-Dirac 
structure, as |1T3| : 



O2 = -l{iuu)idb) + {dd){dh) + {ss){dh)h 

o 

3 _ 

+ -{{du){uh) + {dd){dh) + {ds){sh)}-^, 
8 

- - {du)(ub) + {ds){sb) - {ss){db)}Q 

+ -{3{uu){db) + 3{du){ub) - {ds){sb) - {ss){db) - 2{dd){db)}jg 
8 

= -^-H{3) + ^H{3') - \h{Q) + ^-H{TE) (28) 

The 3, 6 and 15 indicate the SU(3) irreducible representations. The non-zero entries of the 
matrices H(i) in flavor space are ||10|| : 



H{3y = H{3'f = 1 , if(6)}2 = H{Q)f = 1 , H{&)f = H{Q)f = -1 
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H{lh)f = H{lh)f = 3 , H{15)f = -2 , H{15)f = H{15)f = -1 . (29) 

Here 1 = u, 2 = d and 3 = s with the upper indices indicating anti-quarks and the lower 
ones indicating quarks. 

For AS = 1 decays, one has 

O2 = -l{iuu){sb) + {dd){sb) + {ss){sb)h 

o 

3 

+ -{{su){ub) + {sd){db) + {ss){sb)}y 

o 

- ^{{uu){sb) - {su){ub) + {sd){db) - idd){sb)}e 

+ l{3{uu){sb) + 3{su){ub) - {ss){sb) - {sd){db) - 2{dd){sb)}jg 
8 

= -li/(3) + ^HiW) - \h{6) + ^-H{T5) (30) 
The non-zero entries are 1 11 Oil: 



His)' = HiSf = 1 , Hi6)\' = H{6)f = 1 , /f(6)f = /f(6)f = -1 , 

H{TE)\' = H{TE)f = 3 , H{lE)f = -2 , H{lE)f = H(lb)f = -1 . (31) 

For AS = 0, the operators Oi^2, O^^q, and Oy^io can be decomposed as 



where 



3.. 1,. 1^. . 1 
4 

03 = 0-,, 04 = 03' , 

O9 = ^Oi - , Oio = ^O, - hj,. (32) 



O3 = {uu)v-A{db)v~A + {dd)v~A{db)v~A + {ss)v-A{db)v-A 
O31 = {du)v~A{ub)v~A + {dd)v-A{db)v-A + {ds)v~A{sb)v-A 
Oq = {uu)v~A{db)v^A - {du)v-A{ub)v-A 
+ {ds)v~A{sb)v-A - {ss)v-A{db)v~A , 
Ol5 = 3{uu)v~A{db)v-A + S{du)v-A{ub)v-A - {ds)v-A{sb)v-A 

- {ss)v-A{db)v-A - 2{dd)v~A{db)v-A- (33) 
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The operators O5 and Oq have same SU(3) structure as O3 and O4 but different Lorentz- 
Dirac structures. O7, Og have the same SU(3) structure as O9, Oio, but again have different 
Lorentz-Dirac structures. Similarly one can obtain the decomposition of the operators for 
/S.S = 1 case. 

Since we are only concerned with flavor structure in SU(3), operators with different 
Lorentz-Dirac structures and different color structures can be grouped together according 
to their flavor SU(3) representations without affect the results. As long as flavor structure 
is concerned, the effective Hamiltonian contains only 3, 6 and 15. These properties enable 
us to write the decay amplitudes for B PP in only a few SU(3) invariant amplitudes. 



C. SU(3) Decay Amplitudes for B ^ PP Decays 



We will use B, = B^, B^) = (5", 5°, 5°) to indicate the SU(3) triplet for the three 
i?-mesons, and M to indicate the pseudoscalar octet M which contains one of the P in the 
final state with 

\ 



V2^ V6 



M 



TT ' 



K- 



n 

K- 



'V2^ Ve 



9 Vs 



One can write the T amplitude for B PP as 



(34) 



T = AlBiH{'i)\MfMl) + C^BiMiM^H{3y 
+ AlB,H{QyiM]Mf + CjBMjHiey'Ml 
+ A^B^HiTEy^M^Mt + C^B,M;H{TEy'Ml , 



(35) 



due to the anti-symmetric nature in exchanging the upper two indices of H]^{Q), and the 
symmetric structure of the two mesons in the final states, Cq — Aq always appear together 
||T0| . We will just use Cq to indicate this combination. There are 5 complex independent 
SU(3) invariant amplitudes. The results for each individual B decay mode are shown in Table 
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TABLE II. SU(3) decay amplitudes for B ^ PP decays. 



A5 = 

tX- id) = 2Af + + Ci + A^+ 3C§, 
'^R°K° id) = SAg" + Cj - - 3y4Z_ - C^, 

TX(d) = ^,i-cl + ci + 5Ai^+cg, 

Tn%sid) - ^(2Af + \Cl - CJ - Al^+Cg, 
T^'^.-id)=Ci + Ci-A^ + 3C^, 
T^o^o{d)^~^^{Ci + Ci-Al^-5Cl^), 
^K«Jd)^-^^iCi + Ci-Al^-5Cl^), 



AS = -1 

tX- i^) = ^M- + 3^?5 + ' 

T^;^_ (.) ^ + + Ci + AI^ + iCl^, 
Tf/ofio (*) = 2^1" + - - 3Alg. - C-^g., 
^.\(^) = ;^(C^J + 2AI,-2C§), 
T,t.s(^) = V2(4 + M-^§-2C§). 



II. Similarly one can write down the expressions for the penguin induced decay amplitudes 
P. 

Since there are both tree and penguin amplitudes Cj , Aj and Cf, A[, there in gen- 
eral, 10 complex hadronic parameters (20 real parameters). However simplications can be 
made by noticing that Cy^g are very small compared with other Wilson Coefficients, their 
contributions can be neglected to a very good precision. In that case, from Eq. (0), we 
obtain 



"lO r~<Tt- aT\ 3 Cg <-io^T/ 



2 ci - C2 - 3(cr - cl'§)/2 2 ci - C2 

We have checked that the approximation signs in the above are good to 10~^. 

At leading order QCD correction, the above relations are renormalization scale inde 
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pendent, and therefore to this order the coefficients Ci and Ai are also so. This can be 
seen from the fact that when keeping terms which mix only between Oi{Og) and 02(0io), 
the dominant QCD correction gives: Ci(9)(yu) + C2(io)(Ai) = rj^^^[ci{^){'mw) + C2(iQ){mw)] 
and ci(9)(/i) - C2(io)(/.) = ?7"^'''^[ci{9)(mvy) - C2{iQ){mw)]- Here ci,2,9,io(?«vy) are the initial 
values for the WC's at the W mass scale with Ci{iQ){rnw) = 0, r] = as(mvy)/as(/u) and 
= 11 — 2//3 (f is the number of quark flavors with mass smaller than /i). These relations 
lead to (c9(/i) ± cio(/i))/(ci(/i) ±02(^1)) = ±cg{mw) / C2{mw) independent of /x. Mixings 
with other operators and higher order corrections introduce dependence on renormalization 
schemes. We have checked with different renormalization schemes and find that numerically 
the changes are less than 15% for different schemes. Although the changes are not sizable, 
there are scheme dependence. The total decay amplitudes are not renormalization scheme 
dependent, therefore the hadronic matrix elements determined depend on the renormaliza- 
tion scheme used to determine the ratios, {cg ± cio)/(ci ± C2). One should consistently use 
the same scheme. 

Using relations in Eq.(P^, one finds that there are less independent parameters which 
we choose them to be, C^'^{A^'^), Cj, and Cj^{Aj^). Using the fact that an overall phase 
can be removed without loss of generality, we will set to be real. There are in fact only 
13 real independent parameters for B —>■ PP in the SM. 

One can further reduce the parameters with some dynamic considerations. To this 
end we note that the amplitudes correspond to annihilation contributions, as can be 
seen from Eq.(^) where Bi mesons are contracted with one of the index in H{j), are 
small compared with the amplitudes Ci from model calculations and are often neglected in 
factorization calculations P,p6|. Neglecting all annihilation contributions, we then have just 
7 independent hadronic parameters in the amplitudes 

Cf, Cfe''-^ Cje''\ Cy'T,, (37) 

The phases in the above are defined in such a way that all C]"'^ are real positive numbers. 
We will make the assumption that annihilation amplitudes are negligiblly small in our 
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later analysis and leave the verification of this assumption for future experimental data. We 
point out that this assumption can be tested using K~K^, Bg —>■ tt+tt^, tt^tt^ in 

B —>■ PP decays, because these decays have only annihilation contributions as can be seen 
from Table 2 |2|,|T3l. 

D. Constraint on 7 from B —>■ PP Decays 

We are now ready to carry out a analysis using data from B —>■ and B —>■ Kn. 
The experimental data to be used are shown in Table III. 

In general the errors for the experimental data in Table III are correlated. Due to the 
lack of knowledge of the error correlation from experiments, in our analysis, for simplicity, 
we take them to be uncorrelated and assume the errors obey Gaussian distribution taking 
the larger one between ct+ and cr_. to be on the conservative side. When combining from 
different measurements, we take the weighted average. For the data which only presented 
as upper bounds, we assume them to obey Gaussian distribution and taking the error a 
accordingly. 

The analysis in this case is to minimize the given in the below 

where the summation on i is for the available decay branching ratios and CP asymmetries 
listed in Table III. <7Br,cp are the corresponding errors. Here iKib/Kfel) is the due 

to uncertainties in A and \Vub/Vcb\ as in section II. The branching ratios Br{i) and CP 
asymmetries Acp{i) expressed in terms of decay amplitude A{i) = {GF/\/2)[VubV*gT{i) + 
VtbV^*gP{i)] for a particular B P1P2 are given by 

Br{^) = ^ {\A{{)\' + m\')Xp,p, , 

ilTiniBi- B 

^^'^'^ = iA(.)p + |A(.)r 

where Xij = [(1 — (mj + nijY' /m\){l — (nii — /m\)Y/'^ . The amplitudes T{i) and P{i) 
for each individual decay can be read off from Table II. 
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TABLE III. Rare hadronic charmless B — > tttt and B Kn data. The branching ratios are in 
unit of 10"^. 



Br and Aqp 



Data 



Value used from combined data 



Br{B TT+TT-) 



Br{B vr-vrO) 
Br{B K+TT- 



Br{B K-tt"^) 



Br{B K"Tr 



Br{B K°7r°) 



Br{B K-K^) 

Br{B ttV) 
Acp[B K-TT^) 

Acp[B K+TT-) 

Acp{B K^TT-) 



4.3t\i±0.5 [0j28| 
5.9t|f ±0.5 |9| 



4.1 ± 1.0 ±0.7 I 
|2 



J -|^+3. 6+0.9 



3.0-1.2 E^l 
+2.5 



17.2+11 ± 1.2 0,|8| 
16.7± 1.6l|:? 13 
17.0l|o't|° @ 

18.2^^:^ ± 1.6 0,11 

i3.i+ii ± 2.6 m 



14.6111 ±2.7 i 

< 5.1(90%C.L.) 

< 5.0(90%C.L.) I 

9 1+1.7+0.7 
^•-^-1.3-0.6 

-0.29 ±0.23 13 
0.0191°:?^? i2 
-0.04 ±0.16 |3i 
0.043 ±0.175 ±0.021 



0.18 ±0.24 m 



4.4 ±0.9 



6.2 ±2.4 



17.3 ± 1.6 



13.7 ±2.6 

16.2 ±3.8 

14.6 ±4.6 

0.6 ± 1.9 

2.1 ± 1.8 
-0.13 ±0.16 

-0.003 ±0.12 

0.18 ±0.24 
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We use Vcb = 0.0402 ± 0.0019 and iKb/Kbl = 0.090 ± 0.025 in the fitting 
with exact SU(3) symmetry are shown in Figures 5 and 6 by the sohd curves, 
values for the hadronic parameters are 

Cf = 0.13, C^^OM, Cj = 0.13, C§ = 0.16, 
5-3^-27°, 4 = -20°, %=35°. 

^1 

0.9 
0.8 
0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 


-1 -0.8 -0.6 -0.4 -0.2 -0 0.2 0.4 0.6 0.8 1 

FIG. 5. The constraints for p and r) using data from 1 146/146 1 and rare 5 — tttt and B — Kir. 
For the fit with exact SU(3), the best fit value is indicated by the "+" symbol and the X^~Xmin — ^ 
(39% C.L.) allowed regions are inside the region in the solid curve. For the case with SU(3) breaking 
effects, the best fit value is indicated by a diamond shaped symbol and the 39% C.L. allowed region 
is inside the dashed curve. 

And the best fit values for p, rj and 7 are 

p = 0.02, r; = 0.40, 7 = 87°. (41) 

The constraint is weak. We have not given the 68% allowed ranges because to that level, 
the constraints arc basically given by |14()/K6|- We have to wait more accurate data to 
obtain more restrictive constraints. 

At present the errors on the asymmetries are too large and do not really provide stringent 
constraints. However, we include them here hoping that they will be measured soon. By 
then one can easily include them in the fit to obtain more stringent constraint on 7. 
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The results 
The best fit 



(40) 




.^5.8 
5.6 
5.4 
5.2 
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4.4 
4.2 

4 
3 8 

' 20 40 60 80 100 120 140 160 180 

y 

FIG. 6. as a function of 7 using data from l^^b/V^bl and rare B -kit and B Kn. 

The curve (a) is for the case with exact SU(3), and the curve (b) is for the one with SU(3) 
breaking effects. The curves (al) and (bl) are for the cases with the additional condition 

C^e'^s - C| e'*« - = with exact SU(3) and with SU(3) breaking effects, respectively. 

In Figure 5, we show the regions allowed by ~ Xmin = 1 in the p — t] plane by the 
solid curve. As mentioned that at present the constraint is weak which can also be seen 
from Figure 6 where the minimal as a function of 7 is shown by the curve (a) for the case 
with exact SU(3) symmetry, although the Xmin P^^ degree of freedom is smaller than 1. The 
68% allowed region is actually the same as that from iKib/Kftl alone. However, when more 
precision data for rare charmless B — > PP become available, the restriction will become 
more stringent. For example, if the error bars for all the quantities are reduced by a factor 
of 2.45, then the regions in Figure 5 correspond to 95% C.L. allowed regions. 

SU(3) may not be an exact symmetry for B — > PP. We now estimate SU(3) breaking 
effects. The amphtudes for B ^ nn and B Kn will be different if SU(3) is broken. At 
present it is not possible to calculate the breaking effects. To have some idea about the size 
of the SU(3) breaking effects, we work with the factorization estimate. To leading order the 
relation between the amplitudes for B ^ tttt decays Ci^nn) and the amplitudes for B Ktt 
decays Ci^Ki:) can be parameterized as Ci^Kn) = rCj(7r7r), and r is approximately given 
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by 

r ^ ^ = 1.22. (42) 

Here we have assumed that the SU(3) breaking effects in /j and F^'^'^ are similar in mag- 
nitudes, that is, fx/ f-K ~ Ff^~'^ / Fq~*'^ Using the above to represent SU(3) brealcing 



effect, we can obtain another set of fitting results. They are shown in Figures 5 (dashed 
curve) and 6 (curve (b)). The best fit values for the amplitudes are 

C|' = 0.11, C^ = 0.33, Cj = 0.22, Cj = 0.18, 

^3 = 57°, 5e = 200°, % = 85°. (43) 

The best fit values for p, r] and 7 are given by 

p = -0.39, ri = 0.07, 7 = 170°. (44) 

In both exact and broken SU(3) cases, there are two local minimal in the vs. 7 
diagrams. The corresponding values of 7 are very different with one of them around 87° 
and another 170°. These best fit values are dramatically different that those obtained in 
section II. However the best fit values here can not be taken too seriously because, as can 
be seen from Fig. 5, that at 39% C.L. level, almost all allowed range by iV^b/Kfel is allowed 
by i? ^ PP data. At 68% C.L. level, all allowed region by is allowed by data from 

B PP decays. Inconsistence between 7 obtained in Section II and this section can not 
be established. We have to wait more precise data on B PP to decide. 

In the literature it has often been quoted that O12 do not contribute to B^ K^-k' 
and therefore Br{B~ — > K^ti^) = Br{B~^ K^tt'^). In the SU(3) language used here, 
this implies C = Cj e^^^ — Cj^e^^^ — C^e'^^ = 0. This result has been used to derive 
several methods to determine the phase 7. We stress that this is not a result from SU(3) 
consideration and need to be checked. For this reason we also carried out analyses with the 
condition C = 0. For this case, the minimal ^ function of 7 are also shown in Figure 6 
(curves (al) and (bl)) with exact SU(3) and with SU(3) breaking effects, respectively. The 
best fit values with exact SU(3) symmetry are 
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p = 0.05, 77 = 0.41, 7 = 83°, 

Cf = 0.13, C^ = 0.26, Cj = 0.17, C§ = 0.16, 

53 = -13°, 56 = -49°, 5i5 = 27°. (45) 

And the best fit values with SU(3) breaking effects are 

p = 0.12, 77 = 0.39, 7 = 73°, 

Cf = 0.11, C|^ = 0.24, Cj = 0.15, C§ = 0.16, 

53 = -15°, 56 = -56°, 5i5 = 23°. (46) 

The imposition of C = does not force these coefficients to be real. In order to get C 
to be zero, the real and imaginary parts both have to cancel to satisfy the condition. The 
implications of this analysis will be discussed later. 



IV. COMBINED FIT 

In this section we carry out a combined fit of the sections II and III. The total {total) 
is the sum of the x^(II) with sin 2/3 data included from section II plus the x^(III). Here 
X^(III) is the in Eq. (36) of section III with ^^{A, iKfc/Kbl) subtracted. This is because 
that x^(II) already included information from A and iVub/V^cbl- The results are shown in 
Figures 7 and 8. Since x^ill) has a sharper dependence on 7 compared with x^(III), the 
best fit values and errors are dominantly determined by constraints in section II. 

The best combined fit values with exact SU(3) symmetry for the hadronic parameters 

are 

Cf = 0.13, C^ = 0.29, Cj = 0.16, = 0.20, 

53 = -42°, 56 = -20°, 5x5 = 35°. (47) 

In the above we have not given errors for the hadronic parameters because the constraints 
on them are weak. 
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FIG. 7. Constraints on p and rj using combined data from |e_R-|, Ams^^, |14fe/V"cfe|, sin2/? and 
B — > vrvr and B — > i^vr. The three regions from smaller to larger corresponds to — Xmin ~ ^ 
allowed region which is at the 39% C.L., the 68% C.L. allowed region and the 95% C.L. allowed 
region, respectively. The figure on the left is for the case with exact SU(3) and the one on the right 
is for the case with SU(3) breaking effects. The dotted curves are for fit in section II and the solid 
curves are for the combined fit. 

The best fit values for p, r] and 7 and their 68% C.L. errors are given by 

P = 0.I21:IZ V = 0.29tlZ 7 = 67^}^:, (48) 
and the 95% C.L. allowed ranges for p, r] and 7 are 
0.01 < p< 0.29, 0.21 <r] < 0.40, 

-0.86 < sin 2a < 0.45, 0.45 < sin 2/5 < 0.79, 43° < 7 < 87°. (49) 
The best combined fit values with SU(3) breaking effects for the hadronic parameters 

are 

Cf = 0.11, C{ = 0.29, = 0.16, = 0.20, 

^3 = -33°, ^6 = -40°, % = 31°. (50) 

P = 0-inZ V = 0.2n'ol 7 = 661}°:, (51) 

and the 95% C.L. allowed ranges for p, rj and 7 are 
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FIG. 8. as a function of 7 using combined data from \&k\, ^"^b^s' iKtfe/^bl) sin2/? and 
rare B PP data. The dotted and solid curves are for the fit with exact SU(3) and with SU(3) 
breaking effects, respectively. The dashed curve is the same as that from section II with sin 2/3 
included. 

0.01 <p< 0.30, 0.21 K-qK 0.41, 

-0.87 < sin 2a < 0.44, 0.46 < sin 2/3 < 0.80, 42° < 7 < 87°. (52) 

V. DISCUSSIONS AND CONCLUSIONS 

At present the rare charmless hadronic B decay data have large error bars. The main 
contribution to the for the analysis in sections III and IV come from the branching ratio 
for K^n^. In the cases discussed, this mode alone contribute about 2.5 to the x^- The 

best fit value of the branching ratio is only about half of the experimental central value. We 
suspect that there may be some systematic errors in the measurement of this branching ratio. 
If the present central value persists, it may be an indication of badly broken SU(3) symmetry 
or new physics beyond the SM. It is important to improve the precision of experimental data 
to decide whether new physics is needed. 

Because of the large error bars associated with the B PP data, the ranges determined 
for the related parameters have large error bars. Especially the phase 7 has a large range 
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allowed using the B — PP data alone. However, the fit shows no conflict between the flt 
from the consideration using \eK\, ^"^Bd^' I Kit/Kb I and sin 2/5 data. Future experimental 
data will be able to provide a more accurate determination of the phase 7. 

Before closing we would like to make a few comments about our analysis and some other 
related calculations. Our flrst comment concerns the general SU(3) analysis and factorization 
calculations. 

Assuming factorization and SU(3) symmetry, that is the decay constants for all octet 



pseudoscalars P are equal, and the form factors for i? — > P are also equal, one obtains ||26 



= { - ^ « + <o - < + <R) ]X (53) 

where X = f.„F^^'"{mD{m% — m^) and R = m\/{mb — mq){ms + niq). These amplitudes 
are related to the "tree" contributions, a^^ = — with a2i = C2i + C2i-i/N and a2i-i = 

C2i-1 + C2i/N. 

The penguin amplitudes are given by 



C( = -[(4^ + afR) + A(4^ _ af) + l(afo + afR)]X , 



Cfi = + 4o - 4 + 4R)X . (54) 



where a*^ = a* — a^. When small contributions from ay § terms are neglected, one recovers 
the relations in Eq. (PS). Numerically, we have 



= 0.09, ^1^ = 0.42, Cj = 0.26, C§ = 0.15, 
^3 = -15.7°, 6e = -14.5°, % = -14.5°. (55) 

In the above we have used the convention with to be real. 

The amplitudes are in the same order of magnitude as the best flt values in sections III 
and IV, but the phase can be very different. In the factorization approximation calculation 
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here, phases are only due to short distance interaction, rescattering of quarks. Long distance 
contributions can change these phases. The resuhs of the best fit values for the phases 
indicating that there may be large long distance rescattering effects. 

Our second comments concerns the combination of the SU(3) invariant decay amplitude 
C = Cj 6**^5 — Cje"^^^ — C^e^. It has been usually assumed that in the literature that 
C = 0. This leads to Br{B^ K^n^) = Br{B^ K^n^). This result played a crucial 
role in several methods to constrain and to determine the phase 7, for example, using |^ 
B- K-TT^, K°TT-, TT-TT-, B^ {B') ^ K+'R-, 1^+ 1^- [K^ -r' ) , aud B- R-n^, K^tt~, R-?]. 

We point out that C = is based on factorization calculation neglecting annihilation 



contributions and also penguin contributions |26]. In fact, using factorization calculation 



when penguin contributions are included, C does not equal to zero, but C = Cj (penguin). 

(penguin) can be obtained from Eq.(|53|) (the terms proportional to c""^ in Cj. In the 
factorization framework, we can easily check whether C = is a good approximation. Using 
the result in Eq. ( p3D we find that the \Cj [penguin) /Cj\ is of order 5%. It is therefore 
reasonable to assume the penguin contribution to be small and C ~ 0. 

One should also be aware that when going beyond factorization approximation and in- 
clude rescattering effects C may deviate from zero. It should be tested. The fitting program 
proposed in this paper can be easily used to achieve this goal. From the best fit values in the 
previous sections, we clearly see that C can easily deviate from zero. For example, in the 
case with exact SU(3), the best fit value using rare B decay data C is, C = 0.05 — z0.20 and 
with SU(3) breaking effects C is, C = 0.37 + zO.18 which are the same order of magnitude 
as individual Cf . One needs more data to achieve a better test. Until then, the use of the 
methods based on the above equation have to be treated with caution. 

Our final comments concerns the uncertainties in the present analysis. In this paper we 
have developed a method based on SU(3) fiavor symmetry to determine the CP violating 
phase 7. We find that when annihilation contributions are neglected, there are only seven 
hadronic parameters in the SM related to B ^ PP decays. The annihilation contributions 
are small is based on factorization approximation. If it turns out that they are not small, 
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as some model calculations indicated that the penguin related annihilation contribution 
may be sizeable, one needs to include it into the analysis. However, from Table II one can 
see that does not show up in 5 — > Kt: decays, but only to -B — > tttt which is suppressed 
by small Wilson coefficients. One can also carried out an analysis including A^ into the 
fit when more experimental data become available. Future experimental data with better 
accuracy will provide more information. 

In the estimate of SU(3) breaking effects, we have parameterized the SU(3) breaking 
effects in a simple form with Ci{KT[) = {fx/ f'K)Ci{'KTT). In general the SU(3) breaking 
effects may be more complicated. More systematic study of SU(3) breaking effects are 
needed in order to obtain more accurate determination of the phase 7. But in any rate 
we hope that the method developed here will help to provide useful information about the 
hadronic matrix elements and also the CP violating phase 7. 

In conclusion, in this paper we have developed a method to determine the CP violating 
phase 7 based on the flavor SU(3) symmetry. We flnd that present data can already give 
some constraint on 7 and it is consistent with the constraint obtained by using \eK\-i ^it^b^s-i 
|K^6/K^6| ^-iid sin2/3 data. We also carried out an analysis combining data from e^, Am^^^, 
iKift/Kbl, sin 2/3 and data from rare charmless hadronic B decays. The combined analysis 
gives 7 = 67° for the best fit value and 43° ~ 87° as the 95% C.L. allowed range. Although 
there are uncertainties in the fit program, the method developed in the present paper can 
provide useful information about the hadronic matrix elements for rare charmless hadronic 
B decays and the CP violating phase 7. 
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